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Abstract 

In this work we implement the self-consistent Thomas-Fermi model that also incor- 
porates a local conductivity model to an electron-electron bilayer system, in order 
to describe novel magneto-transport properties such as the Drag Phenomenon. The 
model can successfully account for the poor screening of the potential within the 
incompressible strips and its impact on the interlayer Coulomb interaction. An ex- 
ternally applied current in the active layer results in the tilting of the Landau levels 
and built-up of a Hall potential across the layer, which, in turn, induces a tilted 
potential profile in the passive layer as well. We investigate the effect of the cur- 
rent intensity, temperature, magnetic field, and unequal density of layers on the 
self-consistent density and potential profiles of the bilayer system. 
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1 Introduction 



In a two-dimensional bilayer system of charged particles where the separation 
of the layers are large enough to prevent tunnelling, an external current ap- 
plied in one (active) layer can drag the charges in the adjacent (passive) via the 

Email address: siddiki@fkf.mpg.de (A. Siddiki). 



Preprint submitted to Elsevier Science 



1 February 2008 



inter-layer Coulomb interaction, and thus, induce a current. This phenomenon, 
known as the Coulomb Drag (CD) effect, which was proposed theoretically [1], 
and investigated experimentally [2] . The Coulomb drag is of fundamental in- 
terest, since it can be utilized as a probing tool to investigate the screened 
interaction in mesoscopic transport systems. A quantitative measure of the 
Coulomb drag is the transresistivity p2i = E2/J1 where Ji is the applied 
current in the active layer and E2 is the induced electric field in the passive 
layer. 

The presence of perpendicularly applied magnetic field enriches the physics of 
the CD effect by introducing Landau levels. The electrostatic and transport 
properties of the bilayer system were investigated first within the independent 
electron picture [3] and recently within a self-consistent screening scheme [4,5]. 
The main outcomes of the self-consistent scheme are (i) the formation of in- 
compressible strips in the active layer at integer Landau filling factor values, 
which affects the density and potential profile in the passive layer, and (ii) 
the tilting of the total potential within the active layer in the out-of-linear 
response (OLR) regime, induced by a large external current [6]. 

In this work, we combine these effects to formulate the interaction between 
the active and passive layers within the OLR regime, in order to provide a 

self-consistent description of the CD effect. We investigate the effect of the 
current intensity, temperature and density mismatch between the layers on 
the inter-layer Coulomb interaction in the OLR regime. 



2 Model, Results and Discussion 



The electron bilayer system in the present work consists of two GaAs quantum 
wells (having a width of 2(i ~ 3 /xm) that are sandwiched between two silicon 
doped thick (AlGa)As layers. This structure is grown on top of a GaAs sub- 
strate. The layers are translationally invariant. The electron density of the top 
layer is manipulated by a top gate. The electron layers are symmetric (with 
respect to the growth direction) and {h ~) 10 — 30 nm apart separated by an 
undoped (AlGa)As layer. Such a bilayer system is known to be electronically 
decoupled and can be represented by two different electrochemical potentials. 
We note that the electron tunnelling between the layers is not possible. Both 
of the electron channels are formed in the interval —b<x<b, where — fc| 
sets the depletion length. 

The total electrostatic potential of an electron on the line {x, z) due to a 
line-charge at (a;o, Zq) is given by [5] 
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V{x, z) = Vbg{x, z) + Vh{x, z), 

and 



(1) 



2e^ r^r- 

Vh{x,z)^— dxoK{x,xo,z,zo)nei{xo,zo), (2) 

K Jxi 

where — e is the charge of an electron, R an average background dielectric 
constant, and the kernel K{x,xo, z, zo) solves Poisson's equation under the 
relevant boundary conditions given by 



cos^^(x + a;o) +7^ \ 

sin2iL(^_^„)+^2j i ) 

where the z— dependence is given by 7 = sinh(7r|z — Zo\/Ad). The confining 
(background) potential is obtained by inserting a constant number density 
(no) of the background charges into Eq. (2). The gate can be described by an 
induced charge distribution ng{x), that is residing on the surface. In order to 
obtain a flat (gate) potential profile in the bulk we choose the induced charge 
distribution as ng{x) = n°(l + a{x/d)'^) where determines the strength 
of the gate potential whereas a{= 0.7) gives the slope of the induced charge 
distribution. With such treatment of the gate, the average electron densities 
can be modified while keeping the depletion length constant. The Hartree 
potential can be calculated from Eq.(2) for the top (bottom) layer yielding a 
total potential 

V(x, z) = -nj(x, z) - nj(x, z) - Vg(x, z) 



+ VS{x,z)+Vi{x,z). (4) 
The electron densities are calculated within the TFA: 



K{x,xo,z,zo) = - In 



nYix) = / dED{E)f{[E+V{x,ZT,B)-f^T,B]), (5) 

with D{E) is the Landau density of states (DOS), f{E) = l/[exp{E/kBT) + l] 
is the Fermi function, /i^^ is the chemical potential (being constant in the 
equilibrium state) and zt,b is the position of top (bottom) layer, respectively. 
We employ the Gaussian broadened Landau DOS [6]. The self-consistent nu- 
merical calculations start with obtaining the confining potential created by its 
own donors for each layer and the electron densities and the intralayer and 
interlayer Hartree potentials [7] at zero temperature and zero magnetic field. 
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The energies are scaled by the average pinch-off energy Eq — 27re'^nod/R. The 
lengths are scaled by the screening length Gq (= (i*b/2) expressed in terms of 
effective Bohr radius, a*^ = Rh'^/{me^). The electron density and the electro- 
static potential for finite temperature and magnetic field can be calculated 
self-consistently by the above scheme within the TFA. 

The crucial and distinguishing feature of the approach we implement is that 
we neither rely on a one-electron description nor use any phenomenological 
parameters to account for electron interactions, both for inter- and intra-layer. 
In the next step, we drive a relatively high current from the active (here, top) 
layer and calculate the effect of this on the passive layer. We use a local version 
of the Ohm's law to describe current densities employing a Gaussian broadened 
DOS defining the conductivity model [6,8]. We assume a local thermal equi- 
librium, where local quantities vary slowly on the quantum mechanical length 
scales. In the presence of current, it is well known that the Landau levels are 
tilted due to the Hall potential developed [6], therefore the electrochemical 
potential becomes position dependent and thus must be included within 
the self-consistent scheme. This brings two other self-consistent loops: (i) the 
re-calculation of the IJ't,b{^) layer; (ii) the response of the other layer 

to this modification. In this work, we also include this tilting to investigate its 
effect on the potential (see Fig 1) and the density profile of the passive layer. 

Here, we aim to show the effect of the combination of the inter- and the 
intra-layer Coulomb interactions together with the effect of external current 
driven in at least one of the layers. If one drives a current perpendicular to the 
applied magnetic field a Hall potential develops in the transverse direction to 
the current, i.e. Hall potential, resulting from the Lorentz force. If the system 
is all compressible, the Hall potential is linear in ,t, as one would expect to 
observe in a normal metal. Here we observe a similar behavior for the top layer, 
depicted in Fig.l. On the other hand the response of the second layer is as 
expected in the first glance, since, if the top layer has a tilted electrochemical 
potential in one direction the other should be the opposite. In Fig. 1, we show 
the self-consistently calculated potential profiles of both layers. For zero bias 
the top layer (broken black line) and the bottom layer (dotted cyan line) are 
symmetric in x. However, for a large current intensity {Uh-, see Ref. [6]) at 
the top layer to induce an electrochemical potential of this layer (red short 
dashed-dotted lines) a tilt is induced at the bottom layer with an opposite 
slope (dark blue dotted lines). More interestingly, if a finite current is driven 
from the bottom layer with same current direction of the top layer, the tilting 
observed is slightly suppressed (brown dashed-dotted line). 

In summary, by exploiting the smooth variation of the external potential we 
have explicitly calculated the Coulomb interaction in and between the bilayer 
quantum Hall samples in the presence of a strong current driven in one layer. 
Our preliminary results indicate that, a strong Coulomb coupling induces mod- 
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Fig. 1. The calculated potential profiles as a function of lateral coordinate, for zero 
and finite current applied to the top layer. The gate potential is chosen such that the 
resulting density mismatch is at the order of few percents, namely YgjE^ = 0.05. 
Calculations are done at ksT/Eo = 5.5 x lO'^ for no = 4 x 10^^ cm and interlayer 
spacing of /i = 12 nm. The B field is selected such that fuvd Eq = 5.5x 10~^ assuming 
a fixed depletion length h/d = 0.9 

ifications in the potential landscape of both layers, which clearly points that a 
self-consistent microscopic model is necessary to grasp the underlying physics 
of the Coulomb drag effect observed in the bilayer systems. 
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